Spin current and accumulation generated by spin Hall insulator 
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Spin current and accumulation generated by the electric field in spin Hall insulator (SHI) are 
investigated theoretically in terms of the Keldysh formalism. In contrast to the quantum Hall 
system, there are no massless edge modes in general. The spin current is generated near the contacts 
to the electrodes by the hybridization between the metallic states and the conduction/ valence bands 
of the SHI, but is truncated by the sink and source of the spin. However, one can produce the spin 
current flowing out to the conductors, which is attached to the SHI, and also the spin accumulation 
there due to the leakage charge current which breaks the time-reversal symmetry. 
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Intrinsic spin Hall effect (SHE) is an issue of recent 
intensive interests since it might offer an efficient way 
to inject the spins into the semiconductors without us- 
ing magnets/magnetic field at room temperature with 
low energy cost 0, 01- Especially in the hole doped 
case [jj, |2J , this intrinsic SHE is generated by the Berry 
curvature of the Bloch wavefunction in the momentum 
space. Hence the spin current induced by this topological 
property is dissipationless in nature, similar to the charge 
current in quantum Hall effect (QHE). However the volt- 
age drop is produced by the ohmic charge current, which 
produces the Joule heating. Therefore the crucial ques- 
tion to address is " Is there any system showing the spin 
Hall current without any dissipation analogous to the 
quantum Hall system (QHS)?" Related to this question, 
it has been proposed that even band insulators could pro- 
duce the nonzero SHE, i.e., spin Hall insulator (SHI) [4j. 
The idea there is that when the spin-dependent Berry 
curvature has the opposite sign for the conduction and 
valence bands, the occupied valence bands leads to the fi- 
nite spin Hall effect. This is realized when the band gap is 
opened by the spin-orbit interaction (SOI) as in the case 
of HgTe, a-Sn, and PbTe Q . Also the quantum spin Hall 
effect (QSHE) is proposed for graphene , where the up 
and down spins are decoupled into two quantum Hall sys- 
tems. In the generic cases, however, the (pseudo) spin is 
not conserved in the presence of SOI |fj, |j| and the defini- 
tion of the spin current is also obscure. In refs. 0, 0] , the 
conserved spin current is defined by projecting the spin 
operator into each of the band, but this definition can not 
be applied for the spatially dependent situation in the 
presence of the disorder and/or the edges. In refs. 
the another definition of the spin current is given which 
includes the torque dipole density but the procedure of 
fixing the gauge ambiguity is not clear in the present sit- 
uation. On the other hand, the spin and spin current is 
well-defined in conductors without SOI attached to the 
SHI. Therefore it is worthwhile to pursue the spin cur- 
rent and accumulation in the SHI by taking into account 



the realistic experimental setups with the edges and the 
contact to the electrodes/conductors. Especially it is im- 
portant to clarify the difference between the SHI and the 
QHS including the role of edge modes. Here we summa- 
rize the several important differences between SHI and 
QHS: (i) there is no Chern- number characterizing the 
SHE in contrast to the QHE, (ii) the spin current is not 
conserved due to the SOI while the charge current is al- 
ways conserved, (iii) the effective theory of QHE con- 
tains the U(l) Chern-Simons term, which dictates the 
existence of the massless edge modes to recover the local 
gauge invariance for finite size sample while there is 
no such a machinery in SHE of SHI. We will see below 
how these differences affect the low energy properties of 
SHI. 

We start with the following 4-band model for two- 
dimensional SHI defined on the lattice and take the units 
in which K = c = 1 . 
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where T-matrices are defined in ref. Q ■ Let us first con- 
sider the case of A = 0, where the 4x4 T-matrix space 
will be decoupled to two of the 2x2 matrix space. Cor- 
respondingly, the bands are split into the electron- and 
hole-like bands and each of them is still doubly degener- 
ate. When sgn(M) = sgn(t), all bands have zero Chern- 
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number, although they have non-zero Berry curvature 
distribution in momentum space. Therefore, edge modes 
do not disperse across the bulk gap. (It is noted that the 
gap in the edge mode can be narrower than the bulk gap.) 
In the case of sgn(M) = — sgn(f), each band has non- 
zero Chern-number (±2) and there appear edge modes 
across the bulk gap (= 2|M|) just like the QHS as shown 
in Fig. H ( a ) ■ K is interesting that each state near the 
T-point in Brillouin zone is described as a generalized 
Dirac-fermion with quadratic dispersion ||. The spin 
Hall conductivity is quantized in the clean limit, when 
the Fermi energy is in the bulk gap. This corresponds 
to the decoupled two QHS's with opposite sign of the 
Hall conductance and we shall call this case as the quan- 
tum spin Hall system (QSHS). However, it is noted that 
there is a crucial difference between our QSHS and that 
in ref. 5j. In our QSHS, any spin components arc not 
conserved, although the spin z-component is conserved 
in the QSHS of ref. Therefore, we assume the spin- 
independent chemical potential even in a non-equilibrium 
state. This treatment means the scattering between edge 
modes with opposite chiralities on a same edge. There- 
fore, in contrast to the QHS and the QSHS in ref. ||, 
there appears the voltage drop in the edge modes of our 
QSHS. 

Now we consider the case A ^ 0, which causes the hy- 
bridization between previously separated two 2x2 ma- 
trices. This is usually the case in the three dimensional 
model because kinetic terms contains Ti^-matrices. Cor- 
respondingly the hybridization of the edge modes with 
opposite chiralities occurs causing the real gap opening 
as shown in Fig. ^ (b). When M — 0, the bulk gap is 
estimated as 2v2|A|. 
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FIG. 1: Energy dispersion of the semi-infinite system and 
the bulk spin Hall conductance of the periodic system for 
(a) the QSHS, and (b) the SHI. Each curves in the energy 
dispersion corresponds to different fey-value and the dispersion 
is symmetric between k x and — k x . 

Here, we employ the Keldysh formalism applied to the 
finite size sample attached to the electrodes (loj. Re- 
cently this formalism is applied to the 2-band system 
with Rashba coupling ^lj. The real-space Green's func- 
tions are used to obtain the charge/spin current and spin 



accumulation. The effect of electrodes is incorporated 
by the contact self-energy. In this formalism, the lesser 
Green function G < is the most essential correlation func- 
tion to compute physical observables, e.g. charge and 
spin densities, charge and spin current densities. G < 
is defined by Gf ra][r , al] (t, t') = i{cl, a ,{t')c rr7 {t)) , where 

Cra{t) is the annihilation(creation) operator of the elec- 
tron with quantum index a at the position r and the 
time t. For simplicity, we shall focus on the steady state, 
i.e., Gf „ , .At, t') — Gf , r , n (t — t'), and take the 
Fourier transformation with respect to the time index, 

Le -> Gfr^r'a'ft ~ *') G [ra] [r 'a'] ( E ) ' The expectation 

value of observable O r > r — c\,,O r > r c r is given by 
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where Tr^"' means the trace with respect to the index a. 

In order to obtain G < (E), we solve the following 
equations, G<(E) = G R (E)Y* < (E)G A (E), [E - H - 
Y,R{A)( E y^QR{A) = j where the lesser self-energy, 

G r ^ a \E) the retarded(advanced) Green function, and 
Y, R ( A )(E) the retarded(advanced) self-energy. We em- 
ploy the following approximation for the self-energy to 
incorporate the non-equilibrium nature of the system and 
the effects of electrodes and static disorders. 
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where 'E R ' cont (E) is the contact self-energy, and r r (E) 
the local lifetime due to disorders. The local lifetime 
T r (E) is determined self-consistently by the recursion 
equation, - ^ = jN r (E), where 7 represents the 

strength of disorder and N r (E) = ^yTr (<7) [G R r (E) - 

G A r (E)] is the local density of states per unit cell. This 
corresponds to the self-consistent Born approximation in 
real space. The non-equilibrium nature is incorporated 
through the local distribution function f r (E) as 
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= -23?[£ K ' cont (i;)]. f r (E) is determined 
self-consistently by the recursion equation n r (E) = 
^rTr (CT) [G< r {E)\ = N r {E)f r (E), with the boundary con- 
dition at the contacts f r {E)\ nth conta ct = 1+e $^k-^„) > 
where (i n is the chemical potential at the n-th contact. 
It should be noted that the above approximation au- 
tomatically satisfies the charge conservation except for 
the contacts with electrodes where there may be incom- 
ing or outgoing charge current. The higher- level correc- 
tion over the self-consistent Born approximation is taken 
into account through the recursion equation for f r {E) 
which contains the vertex correction of the external elec- 
tric field. Therefore the Ward-Takahashi identity holds 
in the above approximation. 
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The self-consistent equations are solved numerically 
for the system of finite size L x x L y . The setup is 
shown in Fig. El where the electrodes are attached at 
x = ±L x /2, and we take the disorder strength as 7 = 1. 
For the moment, we impose the open boundary condition 
in the y-direction. The difference of chemical potential 
at the two electrodes are fixed at S/i/L x = 5 x 10~ 4 t/a, 

which is in the linear response regime. Here we take 

s 

the definition of the most naive spin current as J r £, — 
h (S^Jrr' + Jrr'Sfj,), where S is the spin-| matrices, and 
J rr i is the charge current. This spin current is not 
conserved in systems with SOI. Therefore there appears 
the source and sink of the spin current represented by 
V • 5 J s » . Fi gure El shows the bird- view of the local den- 
sity of states at the Fermi energy N(Ep), the spin cur- 
rent 5 J Sl * and its divergence V • 5J S » for (a),(c),(e) the 
QSHS (AT = -0.5t, A = Q, // = 0) and (b),(d),(f)the 
SHI (M = 0, A = — 0.5t, no = 0), respectively, where 
fiQ is the chemical potential in equilibrium. It is clearly 
seen that the spin current is flowing near the electrodes 
at x = ±L x /2 in both cases. Here there is a crucial dif- 
ference between QSHS and SHI. For QSHS, the density 
of states N(Ep) due to the massless edge modes is finite 
around the whole edges of the sample (Fig. Eta)), while 
it is produced only through the hybridization between 
the states of the electrode and the conduction/ valence 
bands of SHI (Fig. Etb)). Correspondingly, the spin 
current is terminated at the two ends of each electrode 
(x = ±L x /2, y = zLLy/A) with the source and sink of the 
spin current in SHI, because the edge modes are gapped 
(Fig. |2fd),(f)). On the other hand, the spin current in 
QSHS is flowing from the edges at x = ±L x /2 to the 
other edges y — ±L y /2, although it is not conserved also 
in this case (Fig. Etc), (e)). 

Figure 13 shows the details of the deviations from the 
equilibrium state, i.e., (a) the spin z-component 5S Z and 
(b) the charge current ^-component 5J X in the x = 
cross-section, and (c) the spin current y-component 8Jy 
in the y = cross-section. It is noted that the electrodes 
are attached all through the edges at x = ±L x /2, i.e., 
\y\ < Ly/2. We take the disorder strength as 7 = 1. 
Then, the inverse lifetime is 1/r ~ 0.5t for the bulk and 
0.65i for the edges of the doped systems (no — t), and OAt 
for the edges of the QSHS. Fixing 5fj,/L x = 5 x 10~ 4 i/a, 
we have also investigated systems with other sizes and 
shapes which are not shown in Fig. E] At least in the 
regime 7 > 0.5, 5S Z and SJ X has almost no system size- 
nor shape-dependence. The ^-dependences of SS Z and 
SJ X are also negligible. SJy * has almost no system size- 
nor shape-dependence in the bulk region but scales with 
5fi near the electrodes. This property of SJy' means 
that the sample size dependence of the spin Hall conduc- 
tance is little. The spin accumulation occurs in every 
doped case and the QSHS. This is because the dissipa- 
tion which breaks the time-reversal symmetry is caused 
by the charge current and the associated particle-hole 
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FIG. 2: Local density of states per unit cell at the Fermi 
energy N(Ef), spin current SJ Sz and the divergence of spin 
current V • 5 J s * for the QSHS and the SHI with L x x L y = 
20a x 40a. The electrodes are attached at ±L x /2 from — L y /A 
to Ly/4. The chemical potential of electrons is /io + <5/x/2 at 
x = —L x /2 and fj,o — 8[i/2 at x — L x /2. As for the QSHS, 
the charge current 5 J (not shown) flows along the edges. 



excitations in bulk in the doped case, and at the mass- 
less edge modes in the QSHS. It is forbidden in SHI, and 
hence there appears no spin accumulation, although the 
spin current flows near the contacts with electrodes as 
described above. 

Now we consider the case where conductors without 
SOI are attached at the edges of the sample. Its pur- 
pose is to investigate whether the SHI can inject spins 
into conductors without SOI and make them accumulate 
there. The middle part (\y\ < L y /2) of the heterostruc- 
ture is described by the Hamiltonian eq.(l), and the side 
parts (Ly/2 < \y\ < (L y + L y ond )/2) are the conductors 
represented by the Hamiltonian H — Ylrr' c X^ir'^oCr' , 
where To is the 4x4 unit matrix and t® r , — t for 
r = r' ± ae x<y . Here we take the total system size 
as L x x (L y + L y ond ) = 20a x (20a + 20a). The elec- 
trodes are attached only to the sample part, \y\ < L y /2. 
Fig. 01 shows (a) SS Z and (b) 8J X along the x = cross- 
section, (c) SS Z in the y — —(L y + L y ond )/2 cross-section 
and (d) SJy" along the cross-section a little inside from 
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FIG. 3: Distribution of (a) 8S Z and (b) JJ^ at x = in the 
system with L x x L y — 20a x 40a, and (c) <5J^ Z at y = in 
the system with L x x L y = 40a x 20a. The electrodes are 
attached at ±L x /2 from — L y /2 to L y /2. For clarity, only the 
region of y < is shown in each panel, and SS Z is an odd 
function of y. The curves are shifted in the y-direction for 
5Sz and 8 J x and in the x- direct ion for 5 J x z . 



y = —(Ly + Ly° nd )/2. The inverse lifetime in the side 
conductors is 1/r ~ OM for the doped systems (fio = t) 
and 1/t ~ 0.44i for the QSHS and the SHI, for the dis- 
order strength 7 = 0.5 in all cases. We can clearly see 
the leakage of spins into the side conductors in all cases. 
The discontinuity at the boundary between sample and 
side conductors is prominent. This is attributed to the 
difference of the energy dispersion and density of states 
between middle and side parts. It is noted that spins flow 
out even in the heterostructure of the SHI. The amount of 
the accumulated spins in the side conductors is mainly 
determined by the leakage in charge current. The ob- 
tained spin accumulation is nearly the same as in the 
doped case. Therefore, the SHI is more efficient than the 
doped system, because there is no current in the sam- 
ple in the case of the SHI and there is no energy loss in 
that part. Even in our small system, the energy cost of 
undopcd-SHI is about 45 % of that of doped-SHI. 

In conclusion, we have numerically investigated the 
spin current and accumulation due to the intrinsic spin 
Hall effects in the spin Hall insulator (SHI) by using 
Keldysh formalism. The spin current is flowing near the 
contacts with electrodes, which leads to the spin accu- 
mulation in the attached conductors associated with the 
dissipative leakage current as shown in Fig. [S] 
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FIG. 4: Distribution of (a) SS Z and (b) SJ X at x = 0, (c) 
5S Z and (d) SJ X at y = -(L y + L™ nd )/2,and (e) 5j£* at a 
little inside from y = — (L y + L™ ond )/2 in heterostructures. 
For clarity, only the region of y < is shown in (a) and (b), 
and SS Z is an odd function of y. See also Fig. [5] 
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FIG. 5: Schematic view of the charge current in (a) QHS, and 
the spin current in (b) QSHS, (c) SHI, and (d) heterostruc- 
ture of SHI, respectively. Red and blue boxes represent the 
electrodes with /io + 8fi/2 and /io — <5/i/2 respectively. In the 
QHS, the yellow and light-blue arrows are the deviation from 
equilibrium of electron flow. In the QSHS and the SHI, the 
black arrows are the spin current and the orange and blue el- 
lipses are source and sink of the spin current. The circles with 
dot and cross represent the accumulated up- and down-spins. 
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